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Abstract 

We investigate the evolution of the linear cosmological perturbations in f{R) gravity, an alterna- 
tive to dark energy for explaining the late-time cosmic acceleration. We numerically calculate the 
early-time evolution with an approximation we contrive to solve a problem that commonly appears 
when one solves the full evolution equations. With the approximate evolution equations we can 
fairly assess the effect of the gravity modification on the early-time evolution, thereby examining 
the validity of the general-relativity (GR) approximation that is widely used for the early universe. 
In particular, we compare the CMB photon density perturbation and the matter density pertur- 
bation obtained respectively by our approximation and the conventional GR approximation. We 
find that the effect of the gravity modification at early times in f{R) gravity may not be negligible. 
We conclude that to be self-consistent, in the f{R) theory one should employ the approximation 
presented in this paper instead of that of GR in the treatment of the early-time evolution. 
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I. INTRODUCTION 



The accelerating expansion of the present universe can be explained by an energy source 
of anti-gravity, generally termed dark energy, or alternatively by the large-scale, low-energy 
modification of the gravity theory. In this paper we focus on the f{R) theory of modified 
gravity (for a review, see l|, |2|) with the gravity action 

[ d'x^[R + f{R)]. (1) 



In this theory the deviation from general relativity (GR) is represented by a function of the 
Ricci scalar, f{R), within the gravity action.^ 

The gravitational field equations obtained from the above action are 

(1 + fR)R^.u -liR + f)9,u + {g,uU - V^V.)fR = S-kGnT^, , (2) 

where □ is the d'Alembertian. We use the notation, = df/dR and fjm = d^f/dR^, in 
this paper. 

For the background expansion of the universe we consider a homogeneous and isotropic 
space-time described by the fiat Robertson- Walker metric: 

rfs^ = a2(r) {-c/r^ , (3) 

where a is the scale factor and r is the conformal time. With this metric the above gravi- 
tational field equations lead to 

H^ + H^fn + U-h-\Rfn=^^P. (4) 
6 a 6 3 

(h^_|)+;,(|,^.)_Z_^_£,„.3.G„P, (5) 

where the overhead dot denotes the derivative w.r.t. the conformal time r, the Hubble 
expansion rate H = a /a?, and p and P are the average energy density and the average 
pressure of the universe, for which we will consider matter and radiation. 

For a given expansion history a(r), as well as given p(a) and P{a), Eq. (jl]) becomes a 
second-order differential equation of /(r) or f{R). The functions f{R) that satisfy this 
equation can generate the required expansion history. On the other hand, the dark energy 



^ We consider the metric formalism of f{R) gravity and use the natural units where c — h — 1. 
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models can also generate the required expansion history by choosing an appropriate dark 
energy density Pue{o,) and pressure PdeI^)- Consequently, measurements of cosmic expan- 
sion alone cannot distinguish f{R) gravity from dark energy, and additional independent 
measurements such as the cosmic structures are indispensable. 

For the cosmic structure formation in f{R) gravity, people studied the evolution of the 
cosmological perturbations js-Sj. While the evolution at late times has been widely studied 
{g-S], the evolution at early times is typically treated with a simple approximation, the GR 
approximation, where the deviation from GR is ignored. (For a treatment different from the 

n 

GR approximation for the early times, see [5| where the evolution from the early times to 
the present is studied.) 

In this paper we take into account and carefully investigate the effect of the gravity 
modification in f{R) gravity on the early-time evolution of the linear perturbations. When 
numerically solving the full evolution equations in f{R) gravity for the early times, one 
is usually confronted with a tight-coupling issue. To solve this issue we contrive a better 
approximation, with which we can fairly assess the effect of the gravity modification at early 
times on the evolution. With those at hand, we then examine the validity of the conventional 
GR approximation. In particular, we will compare the density perturbations of the CMB 
photons and matter obtained respectively by our approximation and the conventional GR 
approximation. We will show that the effect of the gravity modification at early times in f{R) 
gravity may not be negligible. Accordingly, for the early-time evolution of the perturbations 
in f{R) gravity, the GR approximation is problematic, and a better treatment is necessary. 



II. COSMOLOGICAL PERTURBATIONS IN f{R) GRAVITY 



For the cosmological perturbations in the early universe we analyze the evolution equa- 
tions of the linear perturbations in the Fourier space and in the synchronous gauge jo]. For 
the metric perturbations we consider the scalar modes, h{k,T) and ri{k,T), defined by the 
line element: 

ds'^ = a^(r) {-dr^ + [5^ + hij{x, r)] dx'dx^} , (6) 
and the Fourier integral: 



kikjh{k,T) + ( kikj - ^^ijj Qv0,r) 



(7) 
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where k = kk and k is the comoving wave number. With regard to the energy part, we 
consider the stress-energy perturbations of cold dark matter (CDM), baryons, photons and 
massless neutrinos. For each of the particle species, 

6{k,T) = 6p{k,T)/p, (8) 
eik,T) = ik^T^,(k,T)/{p + P), (9) 
a{k, r) = -(44 - ls^J)^)ik, r)/{p + P) , (10) 

where S^^ = T'j - 5*jT\/3. 

For organizing the evolution equations of the above perturbed quantities, we introduce 
two new dynamical variables: 

q = h + 6fi, (11) 

X = fRRSRN, (12) 

where 5Rn is defined as the perturbation of the Ricci scalar in the conformal Newtonian 
gauge, and its relation to the metric perturbations in the synchronous gauge is: 

6.. 18H . 4P 1. ^fij 3HH h\ 
6Rn = — V ^^ + — 5 + 3 -r^ + — r^ + — g- (13) 

A. Evolution Equations 

Since the evolution equations of the stress-energy perturbations in f{R) gravity are given 
by the Boltzmann equations with the same form as those in GR [9j],^ here we will simply 
present the evolution equations of the metric perturbations, q, rj, and x- (The information 
about h can be derived from that about q and rj.) 

The i-j component {i 7^ j) of the gravitational field equations in Eq. ([2]) gives 

q = -2aHq + 2P'n - —-j — ^ V p,a, (1 + Wa) • (14) 



^ The Boltzmann equations describe the microscopic physics and therefore the form of the equations is 
independent of the gravity theories. 
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The linear combinations of the 0-0 and 0-j components give 



e (1 + Hf, 
3aH 



2fc^ (1 + fn) ^ 



(. 



H 



Pa9a [l+Wa) 2^ pa^a 



X 



(15) 



X 



aH 



g + /^r/ + 2(l + /^)77+ ai/ - 



1 + fE 



X 



UnGNa'h 
(1 + fn) 



^PaCTa (l+Wa) 



(16) 



In these three evolution equations the subscript a runs over the particle species including 
CDM, baryons, photons and massless neutrinos. In Eqs. (ITSll an d fflGl) the effective energy 
density peS and the effective equation of state Wes are defined as [lO| 



Pes 



VJeS 



SuGn \2 2 a 



PcK 
Pcflf 



hfR " kf + H'^fh 



2a? 



(17) 
(18) 



They characterize the effects of the modification of gravity at the background expansion 
level. 



B. A Problem of Solving the Full Evolution Equations 



One is confronted with a problem when numerically solving the above evolution equations 
for the early universe, particularly Eq. f|T5|) . Here we elucidate the problem. 
We reorganize Eq. ( ITSjl as follows. 



+ 



3aH 



fi 



V 3aH J ^ 3a'^H^ 
^ fR\Q 



V 71— 2^ Pada (1 + Wa) - 



H 

2k* , / k* 
fRV + 



^PaSa 



k^H 



3a?m^ 



3a?m aH^ 



X 



(19) 



This equation can be read as 



f] = [GR-terms + /-terms] / Jr , 



(20) 
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where the "/-terms" denote the terms proportional to the derivatives of / (including x), and 
the "GR-terms" are the other terms that also appear in the GR equations. When the f{R) 
theory is very close to GR (e.g., at early times), and the /-terms are much smaller than 
the GR-terms. In this case, to correctly obtain fj via the above equation, the summation of 
the GR-terms should be as small as //??). However, the error in calculating each GR term in 
the perturbation theory, i.e., due to ignoring high-order perturbations, can be much larger 
than /ij?7, thereby making the calculation of rj in Eq. f|T5|) incorrect. 

This problem is analogous to the issue caused by the tight coupling between photons and 
baryons (before the decoupling around a ~ 10^^) that renders the evolution equations of 
the perturbations in the standard cosmology difficult to solve. A simple but rough solution 
is to invoke the approximation where baryons and photons behave like a single coupled fluid 
with 6.y = 6i). Instead, cosmologists invoke a more accurate approximation, termed "tight- 
coupling approximation", to account for the slip between the photon and baryon fluids jol. 

In our case for f{R) gravity at early times, the largeness of the factor 1//r leads to the 
tight coupling among the GR-terms, which makes the calculation of r] in Eq. (|T5l) incorrect. 
Similarly, a simple but rough solution is to use the GR approximation in Eq. ( lT5l) such that 
the summation of the GR-terms vanishes and the /-terms and furj are ignored. Nevertheless, 
to have better accuracy we contrive an approximation to account for the modification of 
gravity in f{R) gravity. 

The tightness of the coupling can be characterized by aH/fn (for the super-horizon 
modes) or by {k"^ / aH) / f n (for the sub- horizon modes), where the largeness of which suggests 
the tightness. For demonstration we will show in Sec. Illll the evolution of aH/fn in a f{R) 
model, where aH/f^i ~ 10^^ when a ~ 10^'^ and ~ 10^ when a ~ 10~^. 

We know that the CMB observational results are consistent with GR with high precision. 
Therefore the allowed deviation from GR at early times must be small. Accordingly, in the 
viable f{R) models of the late-time cosmic acceleration, fn and its derivatives should be 
tiny at early times, which leads to the tight-couphng issue. 
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C. "Tight- Coupling" Approximation 



For the early-time evolution we construct a new approximation to solve the problem 
discussed above. In dealing with Eq. f lT5]) we decompose rj into two parts, 

r^ = r/(°)+r/W, (21) 

where i]^^'^ and r]^^^ are designed to be comparable respectively to the GR-terms and /-terms 
normalized by H^, so that Eq. f lTSj) can also be divided into two parts which respectively 
lead to the evolution equations of 77*^°^ and r]^^^ with no tight-coupling issue. For this purpose, 
we set 

= p 2^ PaOa (1 + Wa) - Q , f Xp Peff (1 + ^eflf) Q , (22) 

a \ J RJ 

and then derive the evolution equation of r]^^^: 



1 + / 



R 



1 aH ■ aH \ 1 • 



1 / TT Ir \ ^T^GNd fR sr^ s 

-2 «*-TTA;P + p(rT«S''"''"<^+'"°' 



(23) 



Solving Eq. ( 12^ requires the information about x and x- Instead of using Eq. (fTBjl . to 
solve the problem we invoke the following approximation for x'- 

X = IrrSRn ~ X(approx) = SuGNfRRSTN = SnGNfRR ^ Pa(l " 3Wa)5iV,a , (24) 

o 

that is, 6Rn ~ SttGnSTn, where 6Tn is the perturbation of the trace of the stress- 
energy tensor in the conformal Newtonian gauge, 5T/v = {ST^^)n, and 6]\f^a is the density 
perturbation of the a-th fluid in the conformal Newtonian gauge jOj. 

With regard to Xj "we derive its relation to other perturbed quantities from the time 
derivative of Eq. (I2^ 



where 



X ~ X{approx) = StvGn ^ Pa ^/rR " SaH{l + Wa)fRR (1 - SWa)5N,a 

(25) 



4°), ^ (1 + wa) [-Oa-\ + 3^(°) - q 1 , (26) 



That is, we neglect j)^^^ when calculating 5N,a in the above x relation, as the second approx- 
imation. This approximation and that in Eq. (124|) are the two approximations we make in 
our treatment of the early-time evolution. 

To calculate the early-time evolution of the perturbations in f{R) gravity, we solve the 
set of the coupled evolution equations and relations including Eqs. (HM . ( 1211) . ( I22l) . ( 123|1 . 
( ^^ . ( l25l) . and the Boltzmann equations. No tight coupling appears in this set of equations. 



D. GR Approximation vs. Tight-Coupling Approximation 

Conventionally people take the GR approximation jsl Q], where the early universe is 
described by the ACDM model, to solve the early-time evolution equations of the perturba- 
tions in f{R) gravity, thereby giving an initial condition for the late-time evolution equations 
under the matter-domination approximation (and maybe other approximations). In many 
cases people solve the late-time approximate evolution equations from an initial time be- 
tween a = 0.01 and a = 0.03. That is, it is widely believed that the GR approximation is 
valid to a high precision at least before a = 0.01 for most viable f{R) models. 

In the conventional method the effects of the modification of gravity at early times are ne- 
glected and therefore can hardly be assessed. On the contrary, our approximation takes into 
account the effect of the gravity modification in f{R) gravity. Our approximate equations 
in Sec. Ill CI go back to the evolution equations in GR when and fun go to zero, i.e., when 
the effects of the gravity modification are eliminated. Therefore, the GR approximation is 
a limiting case of our approximation and also a rougher approximation than ours. With 
our approximation we can assess the effect of the gravity modification on the early-time 
evolution, thereby examining the validity of the GR approximation. 



III. RESULTS 



We compare the early-time evolution of the cosmological perturbations obtained respec- 
;ively by our approximation and the GR approximation. We modify the CMBFAST code 



ll| to numerically solve our approximate early-time evolution equations of the cosmological 



perturbations in f{R) gravity, while we use CMBFAST to obtain the early-time evolution 
under the GR approximation. 
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FIG. 1: The evolution of —fR/aH, —//Hq , —fu and m = RfRR/{l + Jr) for the designer f{R) 
model with = —1 and the initial condition: fR{ai) = —1.3923 x 10^^^ at = 10^^. 



For the purpose of demonstration, we consider a designer f{R) model [l2j with w^s = 
—1 and the initial condition: fR{ai) = —1.3923 x 10^^^ at Oj = 10~®. This model is 
consistent with the observational results about the cosmic structures 3].^ With regard to 
the other cosmological parameters, we use the values suggested by the Seven- Year Wilkinson 
Microwave Anisotropy Probe (WMAP7) observations [14]: The effective number of neutrino 
species Nes = 4.34, the mass fraction of helium Yhc = 3.26, the Hubble constant Hq = 
73.8km/s/Mpc, the baryon density fraction Q^q = 0.0455, the cold dark matter Q^o = 0.226, 



the effective dark energy Q^m = 0.728, and the matter-radiation equality time z, 



eq 



4828. 



Figure [T] shows the evolution of several /-related quantities for the designer f{R) model 
under consideration, including — f/Hn and the derivatives: —fR/aH (introduced in Sec. 
Ill Bp . —/r and m = R/rr/^I + /r) [1Q|]. The quantity m is conventionally used to charac- 



terize the deviation from GR 



81]. In this model the derivatives of / grow with time from 



tiny values at early times to the order of unity at present, and accordingly / is nearly a 



^ A designer f{R) model with the effective equation of state Weff gives the same expansion history as that 
of a dark energy model with ui^e = Wcs- We invoke the code developed by Wei-Ting Lin to numerically 
calculate f{R) and its derivatives for given Wcs, fR{ai), and the values of other cosmological parameters. 



constant at early times and slightly changes in the recent epoch around the value — 2A. 

We present the evolution of two Fourier modes, k = 0.1 Mpc"^ and k = 0.01 Mpc~^. 
In Fig. |2] we present the CMB photon density perturbation 60 in the synchronous 
gauge and its fractional difference between our approximation and the GR approximation, 
|9o(ours) — 0o(GR)| / [|0o(ours)| + |6o(GR)|]. In Fig. |2]we present the matter density per- 
turbation 6 in the conformal Newtonian gauge and its fractional difference between two 
approximations, |5(ours) — 5(GR)| / [|5(ours)| + |5(GR)|]. The gauge choice for presenting 
6 is made for connecting to the late-time evolution of the matter density perturbation that 
has been widely studied in the conformal Newtonian gauge ^8\. 

In addition, we present in these two figures two relevant quantities: c^, the fractional 
difference between X(approx) and and = {aH/kYm. The fractional difference gives a 
criterion for the validity of our approximation, i.e., the smallness of it indicates the validity 
of the approximation. The quantity Cm is conventionally used to give a criterion for the 
validity of the sub- horizon approximation in f{R) gravity. One may use Cm to determine 
the starting time of invoking the late-time, matter-dominated, sub-horizon approximate 
evolution equations. This starting time will also be the ending time of invoking the GR 
approximation if the initial condition of the late-time evolution is given from solving the 
early-time evolution equations with the GR approximation. While plotting Go(GR) and 
5(GR) from an early time a = 10"^ to a late time a = 0.1, we plot Oo(ours) and (5(ours) till 
the time when c^^ = 0.1 (so as to the fractional difference), before which our approximation 
is valid in assessing the effect of the modification of gravity in f{R) gravity. 

Figure E] shows that for the Fourier mode with A; = 0.1 Mpc~^ the fractional difference in 
the CMB photon density perturbation is about 1% around the photon-baryon decoupling 
time, Zdec = 1090 (a ~ 10"'^), and reaches as large as 10% around a = 10~^. For k = 
0.1 Mpc~^ the fractional difference is about one order of magnitude smaller: < 0.1% around 
the decoupling time; ~ 1% around a = 10^^'^ ~ 0.03. This result indicates that the effect of 
the gravity modification at early times in the f{R) theory may not be negligible compared 
to the accuracy of the CMB observations. With regard to the matter density perturbation 
in Fig. [3l for A; = 0.1 Mpc~^ the fractional difference is about 1% around a = 10"^, which 
is marginally negligible when compared to the current observational accuracy, while for 
k = 0.01 Mpc-i it is smaller: < 10"^ before a = 10~^-^ ~ 0.03. 
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FIG. 2: The comparison of the CMB photon density perturbations obtained respectively by our 
approximation 0o(ours) and the GR approximation 0o(GR). The upper panel is for the case where 
k = 0.1 Mpc~^, and the lower panel for k = 0.01 Mpc^^. We present the evolution of the fractional 
difference in Qo between these two approximations, as well as the fractional difference C;^ between 
X(approx) ^-iid X &s an indicator of the validity of our approximation, and Cm = {aH/k)'^m (where 
m = RfRR/{l + /r)) that is related to the validity of the sub-horizon approximation. 
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FIG. 3: The comparison of the matter density perturbations obtained respectively by our approxi- 
mation (5(ours) and the GR approximation (5(GR). The upper and the lower figures are respectively 
for k = O.lMpc"^ and k = 0.01 Mpc~^. We present the evolution of the fractional difference in 5 
between these two approximations, as well as the evolution of c^^ and Cm- 
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IV. DISCUSSIONS 



In this paper we numerically solve the early-time evolution equations of the linear cosmo- 
logical perturbations in f{R) gravity via an approximation we construct. With our approx- 
imation we can fairly assess the effect of the gravity modification in various f{R) models 
on the early-time evolution of the perturbations, thereby examining the validity of the con- 
ventional GR approximation that neglects the deviation from GR. In particular, we obtain 
the evolution of the density perturbations of the CMB photons and matter, and present the 
factional differences in these two quantities between our approximation and the GR approx- 
imation. This difference indicates the significance of the effect of gravity modification on 
the evolution of the cosmological perturbations. 

We find that the effect of the gravity modification at early times in f{R) gravity may 
not be negligible, particularly for the Fourier modes with shorter wavelengths such as 
k = O.lMpc"^. Thus for self-consistency's sake, the GR approximation is problematic, 
and a better treatment for the early-time evolution is necessary, which our approxima- 
tion may provide. In our demonstration, even though the deviation from GR looks tiny: 
m ^ —Jr ~ 10^^^ when a ~ 10^'^ and ~ 10^^ when a ~ 10~^, the fractional difference in 
the CMB photon density perturbation can reach 1% at the photon-baryon decoupling time 
and even 10% around a = 10^^, which is significant compared to the accuracy of the CMB 
observations. That is, even a tiny deviation from GR at early times may induce a significant 
effect on the cosmological perturbations. This contradicts the conventional thinking. This 
situation is analogous to the issue about the tight coupling between photons and baryons 
before decoupling, which one is confronted with when solving the evolution equations of the 
perturbations in the standard cosmology 9|. 

As a consequence, the CMB observations may provide a stringent test to the currently 
viable f{R) models, meanwhile giving tighter constraints on f{R) gravity than expected, 
and further play an important role in distinguishing f{R) gravity from dark energy. 
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